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The cosmological backreaction arises when one directly averages the Einstein equations to recover 
an effective Robertson- Walker cosmology, rather than assuming a background a priori. While usu- 
ally discussed in the context of dark energy, strictly speaking any cosmological model should be 
recovered from such a procedure. We apply the scalar spatial averaging formalism for the first time 
to linear Robertson- Walker universes containing matter, radiation and dark energy. The formalism 
employed is general and incorporates systems of multiple fluids with ease, allowing us to consider 
quantitatively the universe from deep radiation domination up to the present day in a natural, unified 
manner. Employing modified Boltzmann codes we evaluate numerically the discrepancies between 
the assumed and the averaged behaviour arising from the quadratic terms, finding the largest devi- 
ations for an Einstein-de Sitter universe, increasing rapidly with Hubble rate to a 0.01% effect for 
h — 0.701. For the ACDM concordance model, the backreaction is of the order of Q!^^ ~ 4 x 10"'', 
with those for dark energy models being within a factor of two or three. The impacts at recombina- 
tion are of the order of 10~* and those in deep radiation domination asymptote to a constant value. 
While the effective equations of state of the backreactions in Einstein-de Sitter, concordance and 
quintessence models are generally dust-like, a backreaction with an equation of state Weft < — 1/3 
can be found for strongly phantom models. 

PACS numbers: 04.25.Nx, 95.36.+X, 98.80.-k, 98.80.Jk 
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I. INTRODUCTION 

R has been recognised for some time that the cosmology one recovers from averaging an inhomogeneous distribution 
on a spacehke hypersurface does not in general coincide with a simple Friedmann-LeMaitre- Robertson- Walker (FLRW) 
prescription (see for example [l|, S i, l, i, i, 0, i, i, [3, El, El, Q El, El, E3, El, El HI, HH [H, i23i] for an 
inexhaustive hst). Instead, one recovers correction terms in the averaged Friedmann and Raychaudhuri equations 
which when written as an effective fluid could in principle act as a dark energy. This has two clear advantages over 
the standard picture - it involves no new physics, and the growth of large inhomogeneities expected to cause large 
deviations from the assumed behaviour coincides with the acceleration of the universe. It is with this motivation 
that most studies of averaged and inhomogeneous cosmologies since 1997 have been pursued (see e.g. [24', '25', '2~' 



28, 
58, 



29 
58 



32 
62 



33, 
63, 



34, 
64, 



, M, 39, 4C 

I, m Ei, uL 



42 
72 



43, IM 1451, 146|, |47|, |48|, IM IM [H, ^ [5^, |53|, IM, ^, M 
73j and their references). Recent studies have begun to 
either in perturbation theory (as in, for example, [HI, [gl, [gl, [7l| ) where 
^ for linear modes and grow ing somewhat larger for smaller scales [69| , 
63 , in the latter of which backreactions of the 



quantitatively evaluate the "backreaction' 
the results are consistently of the order of 10 

or in approximate models of inhomogeneous universes, such as in [6 
order of 0.1% are found on scales of lOOMpc"^. 

While it is encouraging that so much attention is being focused on the averaging problem and impact of inhomo- 
geneities in cosmology the cosmological averaging, or "fitting" , problem is distinct from the dark energy problem. 
The prospect of so elegantly solving the coincidence problem is certa inly alluring but the spatial averaging scheme 
popularised by Buchert - or Zalaletdinov's macroscopic gravity [zl, [7l[, [tI, [tI] - is applicable to any cosmology. 
The FLRW cosmology can be obtained by foliating spacetime with maximally-symmetric 3-surfaces. Since in the 
real, inhomogeneous universe these slices are not themselves maximally-symmetric, this involves an implicit averaging 
which should ideally be made explicit. Given a well-defined tensorial averaging procedure, the "cosmological metric" 
would then be that resulting from a direct averaging of the inhomogeneous metric on the 3-surfaces, and may well 
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take the Robertson- Walker form. Since the Einstein tensor is nonhnear, however, the Einstein tensor constructed 
from such an averaged metric wiU not be the same as the averaged Einstein tensor, and the dynamics recovered from 
the average metric and from averaging the Einstein equations will differ. The difference between the two dynamics is 
known as the "cosmological backreaction" . 

Perfect fluid universes (including various scalar-field models) were considered by Buchert [13] but were restricted to 
single-fluid models. An earlier study [lH used a similar formalism but evaluated the effect using multi-fluid transfer 
functions. Studies of late-time backreaction since have tended to focus on pure dust Einstein-de Sitter (EdS) models, 
although 71] included a simplified treatment of radiation alongside the dust components. In a previous study (63j 
(hereafter BBR08), we considered a concordance ACDM universe with ~ 5% baryons and ~ 25% cold dark matter but 
our formalism was not applicable to a universe with radiative or dark energetic components. Backreaction in scalar- 
field universes has also been previously considered by other authors (e. g. psi . [26l. [To] ) but these have focused either 
on inflationary modes or on super-horizon scales. In 2003, Wetterich [245 (henceforth W03) considered a universe 
containing a clumped cosmon field and found a deviation from FLRW behaviour potentially of order unity, with 
an effective equation of state Wcs ~ —1/15; see also [t^] for a study of backreaction on super-horizon scales from 
quintessence with similarly significant results. 

In the present work, we extend our previous formalism to general systems of multiple fluids and scalar fields, 
applicable to realistic inhomogeneous universes. Our formalism allows for the easy incorporation of arbitrary numbers 
of fluids in general metrics and a variety of slicings, and could be applied to a wide range of fully inhomogeneous 
systems in suitable coordinates. For numerical concreteness, we then consider perturbed FLRW models; the formalism 
we have developed makes it straightforward to consider a realistic universe containing some or all of baryonic matter, 
cold dark matter, radiative species, a cosmological constant, and dark energy fluids, across an arbitrary range of 
redshifts. This treatment includes fluid projection terms that were neglected in BBR08. We then evaluate the 
backreaction arising from quadratic terms and its effective equation of state with modified versions of the cmbeasy 
[H, [tI] and CMBFast [zl, [13, [Ml Boltzmann codes, allowing quantitatively concrete evaluations for general models 
without the restrictions of employing, for example, approximate transfer functions as in [T^ . 

Our findings agree with and extend our previous results. The largest deviations arise in EdS universes, with a 
low-Hubble rate model {h — 0.45) generating a backreaction with an effective energy density fi^jj = 5 x 10~^, and 
a model with h = 0.701 an effective energy density of H.'^g = 1.2 x 10^*, comparable to that in radiation. The 
effective equation of state is w^g = 1/57. In deep radiation domination, EdS models with varying Hubble rates tend 
to the same constant impact. The deviations from FLRW behaviour in the concordance ACDM model are suppressed 
compared to the EdS case, with fioff = 4.5 x 10~^ and w^g ~ 1/120. In contrast to an EdS model, in the ACDM case 
the present-day effective equation of state depends on the Hubble rate. The impact at recombination is of the order 
of 10~^, potentially raising the possibility of a detection on the microwave background with a sufficiently accurate 
probe. 

The backreaction from a dark energy model with iti^ = const is similar to ACDM, with a slight increase in fl'^g 
due to dark energy perturbations. The effective equations of state are likewise similar, and increase monotonically 
with decreasing to an asymptote of w'^g « 1/57 for strongly phantom fields. Test cases for dynamical dark 
energy are provided by the standard parameterisation of the equation of state w{a) = wq + {1 — a)wa, an exponential 
potential, an inverse power-law potential and an early dark energy model. The exponential model is a tracking field 
and resembles EdS, but generates a smaller backreaction with Vl^g w 7 x 10~^ and w^g « 1/70. The other models 
are accelerating models and produce a backreaction similar to that of ACDM. Dark energy models generally produce 
a smaller backreaction than ACDM due to the influence of dark energy perturbations, and the effective equations of 
state are always greater than zero. Finally, we consider models with a non-canonical speed of sound, approximating a 
"clumping" dark energy by setting the rest-frame sound speed to zero for a range of constant equations of state w^. 
Doing so, the effective equation of state peaks at w^g w 1/30 for k, —1/2 and declines rapidly for smaller w^. A 
small but accelerating backreaction with vJ^g < —1/3 is found for phantom fields with < —1.87. In combination 
with the results in W03, our study suggests that a more refined calculation is necessary to verify that the average 
behaviour of a universe filled with a scalar field indeed coincides with a universe filled with a homogeneous scalar 
field, and also that sub- horizon backreaction could act as a brake on quintessence models analogous to that in [77| . 

We begin in fJII] with a brief overview of the averaging scheme and write the averaged cosmological equations for a 
general system of multifluids in a relatively general 3+1 split. In i jllll we apply this formalism to a spatially-curved 
Robertson- Walker universe in Newtonian gauge. In ijlVI we first derive some useful approximate solutions in matter 
and radiation domination before numerically calculating the backreaction from quadratic terms in a range of perturbed 
models. We finish with a discussion in ^jV] 
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II. COSMOLOGICAL AVERAGING 



We present only the essentials of our approach and refer readers to [l7|, and BBR08 for further details. Let us 
foliate spacetime with 3-surfaces with vanishing shift vector and the normal to the surfaces 

-,o). (1) 

This system has the line-element 

ds^ = -a'^dt^ + hijdx^dx^. (2) 

An arbitrary matter source with stress-energy tensor T^i, has a density and stress tensor given on the hypersurface 
by 

g = n'^n''T^^, Sij^T^j. (3) 

In a more general scheme capable of averaging tensorial quantities we would include the impacts from the current but 
as this approach considers only the Hamiltonian constraint and the evolution of the extrinsic curvature this is not 
necessary here. 

Denoting an average in a spatial domain V with angle brackets, one may derive [HI, [l^, |T3, from the Hamil- 
tonian constraint and the evolution of the extrinsic curvature an averaged "Friedmann" equation and an averaged 
"Raychaudhuri" equation, 

f ^) ' = ^ {a'^) + (a^) ^-liQv+n^), '^ = -'^{a^{g + S)) + (a^) ^ + \{Q^ + V^) . (4) 
\aT> J 3^'^'36 a-D 3^ '^'33 

These are often known as the Buchert equations, although it should be noted that they have been rederived multiple 
times with different notations (e.g. P, UjI)- The kinematic and dynamic backreactions and the curvature correction 



are 



Qv = {o? (k^ - K]k() ) - ^ {aKf , Vv = (aD'D.a) - (aK) , TZv = (a^U) (5) 
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where Di is the covariant derivative on the 3-surfacc, TZ the Ricci scalar on the 3-surface, and 

1 

2a 



Kij = - — htj (6) 



the extrinsic curvature. The kinematic backreaction is related to the variance of the extrinsic curvature, the dynamic 
backreaction vanishes in a synchronous slicing, and the curvature correction is a direct average of the Ricci 3-scalar. 
The Hubble rate is defined by the volume expansion as 

av o V 3 

From these equations (jj]) one may further derive an integrability condition connecting the backreaction and curvature 
terms. In a general slicing and with general fluids its form is rather complicated (it is presented in BBR08), although 
it simplifies significantly for a synchronous slicing when the dynamical backreaction vanishes, and more so for dust 
models (see [l3| for details) for which it is of significant use [H, . 

We consider a universe filled with baryons, cold dark matter, photons, three species of massless neutrinos, a 
minimally-coupled scalar field, and a cosmological constant. The perfect fluid provides us with a good description of 
all of these. The underlying model treats radiation in the usual manner as effective fluids governed by Boltzmann 
hierarchies and so no physical error is introduced in doing so; the averages are purely formal and we can unambiguously 
map the lower moments of the hierarchies onto the fluid density and velocity. A perfect fluid with 4-velocity 



(8) 



normalised to u^u^ = — 1, rest- frame density p and rest-frame isotropic pressure p has the stress-energy tensor 

Tf^^ = pu^u^ -I- pP^^u, Pf^u = g^ii^ + u^u„. (9) 
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The projection tensor P^j^ projects quantities into the fluid rest-frame. The forms for baryons and CDM are then 
recovered by setting p = 0, those for radiative fluids by setting p — p/i and that for a scalar field with p = w^pp. The 
fluid contributions to the average cosmological equations for a species (a) can then be written 

^T^G , 2 \ SttG / 2 \ I 4:nG , 2 c \ A / 2 \ l-7r(Q) /1n^ 

[a 0{a)) ^ [a p(a)) + J^v ^ ^ (a •S(a)) = -47rG'(a p) - -J^^' (10) 

where 

. ^ (,.,^ _ ^,^^^) . ^ (,.^ _ p.^) T^H) (11) 

accounts for the tilt between the fluid rest-frames and the foliation. 
The multifluid averaged cosmology can then be written as 



a \ a 



(13) 



It is tempting to immediately define an effective energy density, pressure and equation of state of the "backreaction 
fluid" 

8ttG 



-P, 



3 .off = J2^v-liQv+n-o), (14) 

a 

167rGp,ff = ^^Ses = 2j2^v +liT^v~3Qv-4Vv) (15) 



and so 



WeS = -77 7f^ ■ (16) 

(Compare with, for example, equation (45) of BBR08). Note, however, that while this does give a correction to an 
effective FLRW it does not let us directly compare with an input FLRW model: 7?. 15 contains a background curvature 
and, depending on the choice of a, Vv could contain background functions of the Hubble rate. Furthermore, additional 
non-linear contributions arise from the fluid and cosmological constant terms. Accordingly, the effective fluid is defined 
in a model-specific manner. 



III. NEWTONIAN GAUGE PERTURBATION THEORY 



The system of equations presented in ^|TI]is general for any system for which the 3-1-1 split remains valid. To make a 
quantitative evaluation of the impact of the backreaction, we must specify a system. The simplest first approximation 
is a perturbed Robertson- Walker metric, which as in BBR08 we consider in Newtonian gauge. The use of perturbations 
and this gauge is subject to the same caveats as before. Firstly, linear (or mildly non-linear) theory restricts us to 
averages on large scales, which in general will be outwith our past light-cone and not necessarily observable. Secondly, 
using perturbation theory automatically assumes that the backreaction terms are small - otherwise, the background 
would be poorly chosen. While this is undoubtedly so, and while it cannot naively be taken to give us the impact 
observed on the CMB (although |72| attempts to connect a similar formalism with CMB observables) or in supernovae 
surveys for which averages across much smaller scales (such as those in [68l [6^ [83|) must be performed, the use of 
perturbation theory means that we can quantitatively evaluate backreaction terms in the simplest realistic cases, 
which can provide us with some of the features of a more general model. Moreover, linear theory is still valid on 
scales above k ~ 0.06/iMpc~^, as seen from the matter power spectrum [sl], and at the epoch of recombination and 
in radiation domination the universe is extremely well described by perturbation theory. Second-order theory, even 
in the present epoch, is valid to slightly smaller scales. In these cases our calculations will be robust. 

It is also worth commenting on the different scale factors we employ. When we use a perturbed Robertson- Walker 
universe as the underlying model the scale factor a(t) is that associated with the metric, appearing in and governed 
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by the Einstein equations. The scale factor aj}{t) appearing in the averaged cosmological equations is, in contrast, 
a "reconstructed" scale factor, defined by an average of an inhomogeneous distribution across the domain T). Its 
evolution is governed by the averaged equations and clearly has no impact on the underlying model. The distinction 
between the two is extremely important. 

This issue arises because of the somewhat circular nature of the argument: we are assuming a large-scale "average", 
adding perturbations to approximate the inhomogeneous universe and re-averaging the results. Phrased differently, we 
are defining perturbations on a smooth "background" manifold, and then averaging across the inhomogeneous manifold 
defined by the perturbations. In a more realistic model, one would not be employing an underlying Robertson- Walker 
model at all but rather a distribution of sources (as in [68|) or multi-scale models (as in [l^ll^]), and the only scale 
factor would be the reconstructed ax), and this confusion would not arise. 

Take the Newtonian metric in the form 



where the spatial metric is 



ds^ = -{! + 2-^)dt^ + a^(i)(l - ^^Hjdx'-dx^ (17) 



7., = diag ( ^-^,r2,r2sin2(0) 1 (18) 



and we work with spherical polar co-ordinates = {r,9,(j)). The perturbations ^' and $ are non- linear and can if 
required be expanded in a Taylor series; we restrict ourselves to second-order, retaining up to quadratic terms in 
perturbations. Therefore 

* = E^*(»)^*(l) + ^*(2)' * = E;^*(")=*(l) + ^*(2)' $* = <i>(l)*(l) (19) 

n n 

with similar results for $^ and ^P^. We then immediately have 

a2^1 + 2*, a = 1 -I-* - i^-^ + 0(4'^), a"^ = l-* + ^*2 + C'(*^), ^ 1 - 2* + 4^'^ + C'(5'3) (20) 

and 

= a\t){l - 2$)7„-, h'^ = a-\t){l + 2$ + m^)^' + 0($'). (21) 
Following BBR08, the extrinsic curvature is 

= -(--{1 + 2$)$) 6} (22) 



and so the kinematical backreaction is 



Qv = 6{{^'}-{^) ), (23) 

trivially related to the variance of $. Expanding the covariant derivative on the Newtonian slicing and using the time 
derivative of the lapse, we find the dynamical backreaction 

Vv = ^ (V^* - (V«')2 + 2$V^* - (V$) • (V*)) + 3- {^'i' - 2*^'^ - 3 (i>^^ 

^ (— + 2$—^ + (— + 2$— + — ^ — ^ \ (24) 

dr dr J \ dr^ dr'^ \ dr J dr J / 

where V — {d/dr, {l/r)d/d9, {1 /r sin 9)d/d4>) is the usual Euclidean gradient operator and the curvature corrections 
are written explicitly. The Ricci scalar on the 3-surfaces generates the curvature term 

7^p = ^ + ^ (2V2ci> + + 4(2$ + v|/)v2$) + 1^ + - - 

a a \ or or'' 

4(* + 2$)r2^ + 3rM ^ + 12 (* + 2$) - 12$ (* + $) ) . (25) 
a"^ \ ar^ \ Or I or / 
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As expected, this contains the background curvature from the input FLRW model. 
The fluid 4-velocity is 



(26) 



and its norm is 

hijv'v^ = a^V^ = (j^j^ + r-'vl + sin^ Ovf^ . (27) 

We linearising the density and the pressure with respect to the FLRW background, 

p = + p = p(l + 7), (28) 

where again 8 and 7 contain contributions from first- and second-order. This then gives the total fiuid contributions 
from this species as 

— {a^p)+J^Tj = —^ + Tv, —{a\p + S))+Tv = —{-p + ^P) + ^T-D+ST, (29) 

where the corrections to the FLRW case have been consolidated into a density correction and a pressure correction 
S-D, defined by 

O'T' O C 

1^ =p(5 + 2* + a2(l + uJ)V2 + 2*^), ^=^{2.cl8 + &w^ + o?{l + w)V^ + &cl^5) . (30) 

Setting /C = 0, M) = = 0, neglecting anisotropic stresses (implying $ = \1/) and swapping to Cartesian co-ordinates 
reduces these results to those in BBR08. 

The effective fluid that arises in an FLRW universe perturbed to second order is therefore given by 

^75,3. = + 1a (a^ - 1) - i (Q^ + n-D) , (31) 

167rGp,ff = 4<SI,-2(a2_l)A+i(7^I,-3QI,-4p2?) (32) 

where it is understood that TZj) docs not contain the background term. We will characterise the "backreaction" with 
the redundant set of dimensionless quantities, {Qeff, Ai?/it!, Weff} where 



AR _ 4nG (p,ff + 3^,ff) _ ( V3) {Qv + Vv) - Ea (^4"^ + 4"' 



R |a/a| |a/a| 



(34) 



_ p,s _ 1 - 4Pp - 3Qp + 12.$p - 12A (g^ - l) 

fieff and AR/R give the corrections to the Friedman and Raychaudhuri equations respectively that an observer in 
the perturbed FLRW universe would reconstruct. Normalising the correction to the Raychaudhuri equation in this 
manner introduces an artificial singularity if the FLRW model passes from deceleration to acceleration (or vice- versa) ; 
if desired this can be avoided by normalising instead to the contribution from, for example, dust matter. 

IV. RESULTS FROM QUADRATIC MODES 

When averaging across very large scales, due to their Gaussian nature it is a reasonable approximation to neglect 
averages of pure first-order perturbations. The same argument cannot be applied to second-order perturbations. 
However, it is not clear how in general an average of the form ((5(2)) can be performed, although analytic solutions 
could be employed in particular, highly simplified, cases. Here we consider solely the impact from the quadratic terms 
and leave the general issue to future study. While a calculation from purely quadratic terms is strictly neither a 
complete nor a consistent calculation, it does give a firm order-of-magnitude estimate of the effect. Since second-order 
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quantities can be written using the Einstein equations as combinations of first-order quantities, the error introduced 
is at most approximately of order 0(1) and calculations from the quadratic terms treated loosely as upper bounds on 
the impact. 

If we neglect the averages of pure first-order perturbations and their derivatives, which is a reasonable assumption 
when averaging across very large scales, the backreactions and curvature correction reduce to 



\ {2^VH - (V*)2 - (V$) • (V*)) - 6- - 3 



-— ( 6r$— 

or 



a 

6/C 2 



— 

Qj.2 



9$ a* 

dr dr 



a' 



(3(V$)2 +4(2$ + *)V2$) 



'4U^ + 2^)(,r'-l 



a 



dr 



dr^ 



+ 3r^ 



dr 



12$(* + $) 



and the fluid corrections to 



(36) 



(37) 



(38) 



(39) 



4"' = ^ {a+^ia^Wvl, -f 6c^(.,*,(.)) = Ix^^ + ^ {3c%, - 1) (m,^,) . (40) 



Qv 

(a) 



6 / Vi,(k) 



On very large scales where we can invoke the ergodic principle, we can turn the volume averages into ensemble 
averages. If P^(fc) is the primordial power spectrum of 4", then with /C = the various correction and backreaction 
terms become one-dimensional integrals across the primordial power spectrum: 

2 cifc 
T' 

-3 j r^{k) Q*** + + 1 1*1') + I + 

J kV^{k) (5 |$|^ + 2^-$* + 2**$) dk, 



3 

47rG 



P(a) / (3c^(a)*^(a)+3c^(„)'J'*^(a)+(l+'«(a))a'h(a)|') 



2\ dfc 



(41) 
(42) 
(43) 
(44) 
(45) 



The power spectrum is defined in the usual manner as 



27r2 



(vl/(k)**(k')) = ^7'^(fc)5(k-k') 



and we assume a power law form 



r^ik) = As ( — 



(46) 



(47) 



Unless otherwise mentioned we assume the WMAPV parameters Ag — 2.4 x 10^^ and Ug — 0.96 at a pivot fc, = 
0.002Mpc~^ [s^l. Due to the form of the integrals, modifying the amplitude and spectral index will amount to 
a renormalisation of the backreaction. In principle, the expressions (|4m45p are not accurate; this is because the 
Fourier modes are defined with respect to the FLRW background, while the averages are defined with respect to the 
inhomogeneous background. However, the correction term can be shown to be negligible for the models we consider. 

The curvature terms in equations (|37p and (|38p are non-trivial. Our approach involves computing the perturbations 
as a function of wavenumber k and converting volume averages on large scales into ensemble averages. Doing so, we 
can convert gradients into powers of ik. However, the curvature modifications involve derivatives with respect to a 
radial coordinate r and we cannot deal with this so simply. Rather than pursue our current approach in the curved 
case, we consider only flat models with /C = 0. Employing a somewhat different approach, Rosenthal and Flanagan, 
however, recently considered the backreaction in certain closed models 67]. In forthcoming work we will consider the 
spatially-curved cases in our formalism, employing an alternative direct spatial averaging approach. 
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1. Analytic Approximations 

We can derive analytic solutions to the above equations in particular regimes, the most interesting of which are a 
universe filled with a single radiation fluid, which serves as an approximation for the real universe in deep radiation 
domination, and a universe filled with dust, corresponding to an EdS universe or the ACDM model for z > 1 — 2. In 
both cases the anisotropic stress vanishes and the curvature correction and backreactions become 

nv^-^ f kV^{k)m'dk, Q^^l fv^{k)\'f'f^, -pv^lnv-^Qv--- I V^{k)Re{^^')'^. (48) 

(X J a J K- O Qi (X J rv 

Here ' ~ d/drj and rj is the conformal time drj = dt/a. For a single fluid with v'^ = = ''^^ + ''^^ + ; 

Tv = ^ lv4k){2Re{M*) + {l + w)a\')^, Sv = ^Tj, + ^ {3cl ^ l) J V^,{k)Re{M*)^. (49) 

The evolution of the quantities <&, S and v is given by the Einstein and conservation equations. See for example 
[ssl [131 for these; here we will present only the results. 

For a universe filled with a radiative fiuid where w = = 1/3 the background evolves with respect to conformal 
time as 



and the pressure correction becomes 



a' I a" a ri , ^ 

- = — = 0, — = (50) 
a rj a Ooq rjcq 



Sv = \Tv. (51) 



The perturbations can be seen ([86j) to evolve on sub-horizon scales {krj ^ 1) as 

av{k,r^)^^Dsm{kr]), S{k,Tj) = 2D cos{krj), <^>{k,rj) ^ (52) 

The matter terms are dominant and taking = 1 for simplicity implies 

^ Hln^D^ sin^x, 

Tv - 2Sv = -2-^ / dx (53) 

where x — krj. The effective energy density ([33]), modification to the Raychaudhuri equations (|34p and effective 
equation of state ([^S)) are then 

rJoff ~ a^V^Tv, ^ ~ -^cs, ~ ^- (54) 

il o 

Setting ?7eq from the present day using (|63p . the effective energy density is 

f^off ~ 4:D^^Im J ^ dx ^ const. (55) 

The modifications in deep radiation domination, then, tend to a constant and have no dependence on the Hubble 
rate. They will, however, have a dependence on the scalar spectral index. We can also see from ((55|) that, perhaps 
contrary to expectation, the oscillations present in the pre-recombination plasma and in radiation domination are not 
expected to have a significant impact on the results. 

The opposite regime of interest is a universe filled with dust. Setting w ~ — gives 

S^ = ^{a^vl,)>0, (56) 

which corresponds to the "gravitational pressure" evaluated in W03 and BBR08. This term was otherwise neglected 
in BBR08, a good approximation when considering the impact on the Friedmann and Raychaudhuri equations but 
less so when considering the effective equation of state. In a dust universe, the perturbations evolve as 

$(A,r7) = 4'(fc,r;) = $i(fc), $'(fc,ry) = 0, ,5(fc, r/) = -i$i(fc)fcV, oi;(fc, r;) = i$i (fc)fc?7 (57) 
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with the background 



The correction terms are then 



2 



2 a" _ 2 _ 77' 
— , — ij , d — ^ . 



(58) 



2 18 r 4 1 

Qi, = 0, rv = -nv, TZv = — ^ kVi,{k)<^>l{k)dk, Tv = —nv, Sv = --Tv. (59) 
9 a J 81 4 

These exact results agree with the purely numerical results from BBR08, where TZv/Vv ~ 4.5 and T-p/TZv ~ 1/20. 
The effective equation of state is 

WcS = ^ (60) 

The gravitational pressure thus makes the effective equation of state from linear perturbations greater than zero - 
neglecting St> recovers Wcs = —1/19, as found before. The effective energy density and change to the Raychaudhuri 
equation are 

"-ft ^ ?77?2« / kr^ik)^l{k)dk, ^- (61) 

The backreaction in matter domination thus evolves linearly with the scale factor. Since the impact at the present 
day is expected to be of the order of Oj^^j « 10^'^, the backreaction at last scattering is ^leS ~ 10~®. Our calculations, 
arising from linear theory, can only be interpreted at the present epoch as an order of magnitude estimate of the 
backreaction from large-scale modes, and the issue of backreaction from highly inhomogcneous structure is necessarily 
left unaddressed; at the last scattering surface, however, the universe is well-described by linear perturbation theory 
and our result is robust. The level of deviation from Robertson- Walker behaviour at recombination is then similar in 
magnitude to the standard polarisation anisotropics, suggesting it may be observable with a suitably sensitive probe 
and a method of mapping the spatially-averaged cosmology onto observable quantities. The question of impacts from 
backreaction on the CMB was recently discussed in [t^I using scaling assumptions for the correction terms. 

As in the radiative case, a change in the primordial power spectrum will merely rescale the modifications by a 
constant factor. However, in contrast, they are sensitively dependent on the Hubble rate. While might seem to 
imply that a higher Hubble rate will induce a lower backreaction, this is not the case. <&i(fc) is set at the epoch of 
matter-radiation equality and by ([5^ on small scales is approximately of the form 

^m^-A'-^ (62) 



where 



jyoq ~ aoq»7o ~ (63) 

-"0 



This scales (|6T|) by a factor of Hq. Furthermore, we choose to take as our domain the Hubble volume, and so 
fcmin ~ t^Hq at the present epoch and the limit introduces further factors of Hq. Wc characterise the result by 

f^eV « h^'"^ (64) 

where Hq = h x lOOkms^^Mpc^^ and we expect m to be small. The dependence on the Hubble rate is then not 
necessarily trivial, but we expect an increase in the Hubble rate to increase the signature from backreaction. 



2. Einstein- de Sitter and ACDM 



In BBR08 we considered numerically an EdS and a ACDM universe for z < 100, finding Wcs ~ —1/19 at the 
present day, declining slightly as z ^ 100. However, we neglected the "gravitational pressure" term St> and were 
constrained to low redshifts and pure dust models. Here we update the treatment to the WMAPV concordance 
cosmology J^, Isl], high redshifts and dark energy cosmologies. We use modified versions of the cmbeasy [tI] and 
CMBFast [79[ Boltzmann codes (collectively referred to as Backfast) and numerically integrate equations (|4m45p as 
a function of time. Our domain is fc S (27r/7y, 40Mpc^^) where the small-scale limit is set such that the integrals 
converge. Naturally we do not claim that perturbation theory actually applies on such scales, merely that our results 
give the total contribution from such modes. 
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FIG. 1: Effective energy density of backreaction and modification to the Raycliaudhuri equation for EdS and ACDM. Dashed 
lines are negative. Left: QcS for ACDM with Us = 0.96 and Us = 1.04, and an EdS model with h — 0.45 and h = 0.701. Right: 
— AR/R for the same models at high redshift. The impacts are in agreement with the approximate predictions. 



Figure [T] shows flcff and \AR/R\, clearly showing the linear dependence at low— z. The net deviations from FLRW 
behaviour at the present day in a sample EdS universe with h = 0.45 are ft'^g w 4 x 10~^ and \AR/R\ w 3 x 10^*^, 
while those in the WMAPV concordance model are f}°fj « 4.4 x 10"^ and \AR/R\ « 4 x 10"^. For h = 0.701, an EdS 
universe shows a much greater impact with an effective energy density « 1.2 x 10^*^, larger than that in radiation, 
« 4.2 X 10~^/i~^ « 9 X 10~^. Viewed differently, this is a 0.012% effect. Increasing the scalar spectral index from 
Us = 0.96 to Us = 1.04 increases the effective energy density for ACDM to QcS ~ 5.6 x 10^^, with \AR/R\ « 5 x 10^^. 
The ratios between the corrections agree well with the analytic approximations and in all cases the corrections serve 
to decelerate the expansion. 



Figure [5] shows the present day effective energy density fi^jj for an EdS model with a range of Hubble rates and 
all other parameters held constant. It should be noted that these models are illustrative only, and that without 
adjustments to the primordial power spectrum and matter abundances they cannot fit the CMB, let alone large-scale 
structure. The power law 

f7^ff cx h^-^^^^ ^ m = 0.1875 (65) 

fits our results for the range h S (0.4,0.8), which covers most values of interest. 

As expected, the impact at the epoch of recombination is of the order of 10~* and in EdS cases is larger still. Starting 
soon after recombination the corrections interpolate smoothly between the matter and radiation approximations 
presented above. The delay between recombination and the EdS behaviour is due to the photon velocities which 
are decaying but remain significant for a period after recombination. As predicted, the corrections tend to an 
approximate constant value deep in radiation domination, relatively independent of the Hubble rate. If one considers 
the z > 150, 000 regime there is a further turnover and the backreaction begins once more to decline. 



Figure [3] shows the corresponding effective equation of state for the EdS and ACDM models (see Figure [7] of BBR08 
for comparison with the previous results). For the EdS models, w^g = 1/57 as expected. For ACDM the transition at 
z w 0.7 is clearly visible, and the equation of state declines significantly as one approaches the current epoch, ending 
at Wgfj ~ 1/120 for the concordance model. At high redshifts, the effective equation of state rises towards one third 
and asymptotes at WcS ~ 0.31 in both EdS and ACDM models. 
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FIG. 2: Dependence of the backreaction on the Hubble rate for EdS models with Us — 0.96. Left: Q'^g as a function of Hubble 
rate. Points show the numerical results, while the curve is QcS tx h^'^^^"' . Right: The convergence of the modification to the 
Friedmann equation to a constant value deep in radiation domination. 
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— ACDM, h=0.45 

— ACDM, h=0.70 

— EdS, h=0.45 

— EdS, h=0.70 
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FIG. 3: Effective equation of state of the modifications in EdS and ACDM cosmologies. The asymptotic values of WcS — 1/57 
in matter domination and Wcg « 1/3 are clearly visible, as is a transition at z ~ 1000 from the decay of radiative velocity. 



3. Dark Energy and Quintessence 



The main dynamical alternatives to a standard ACDM model are quintessence cosmologies in which a scalar field, 
minimally-coupled to gravity, drives an acceleration of the universe at low redshift. In contrast to the ACDM model 
dynamical dark energy models contain dark energy perturbations, which contribute non-trivially to Tp and iSp, and 
alterations to As we are working on linear scales we still expect ^ a-nd the dominant differences from 
ACDM to come from the dark energy perturbations. Since T-p and S-p contain terms of the form V{{(l))) — {V{(j>)) 
one might expect relatively large deviations, in models with an exponential potential for example. Indeed, in W03 
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FIG. 4: Left: Effective energy density in backreaction and modification to the Raychaudhuri equations in ACDM, a quintessence 
model with Wif, — —0.9 and a phantom model with w — —1.1 (solid curves). Also shown (dotted curves) are are models where 
the rest-frame speed of sound squared of the dark energy is set to zero. Right: The effective equations of state for these models. 

deviations on the order of unity were estimated from a nonlinear, clumping model of such a field. We consider the 
following dark energy models: 

• a constant equation of state, = const, 

• the standard examples for scalar field potentials: an exponential and an inverse power law, 

• two parameterisations of the dark energy evolution, first the standard parameterisation of the equation of state, 
w{a) = wq + Wa{l — a) , and a parameterisation of the energy density which includes an early dark energy 
component. 

Finally, we will investigate the effect of a change of the sound speed of the dark energy component.^ Unless 
otherwise noted, in this section we employ WMAP 5-year parameters h = 0.701, Ug = 0.96 and Q,^ = 0.721. 



Dark energy models with a constant equation of state models serve as a good first example and are indistinguishable 
from ACDM for z > 10. We separate models into the "quintessence" regime with > —I and the phantom regime 
with Wcf, < —1, in both cases taking = 1. In Figure |3] we plot the modifications and effective equations of state for 
models with Wtj, = —0.9 and = — 1.1. In both cases dark energy perturbations make the modifications larger than 
those in ACDM, with the increase in St> making the change in AR/R larger than that in Jleff- The quintessence field 
with = —0.9 in particular generates a large deviation from the standard Raychaudhuri equation. The effective 
equations of state in both cases are similar to ACDM, which occupies an envelope between the quintessence and 
phantom curves. For = —0.9, WcS is slightly less than that for ACDM, while for = —1.1 it is slightly larger, 
and ACDM occupies an envelope between the two curves. We then expect the effective equations of state from models 
with Wif, > —1 to be slightly less than those of ACDM and those from phantom fields to be greater. 



In this study, unless noted we take = 1 for a scalar field cosmology. However, the use of this in our context has recently been 
challenged 18911 . 




FIG. 5: Dark energy models with a constant equation of state and sound speed — 1 (blue) and — (red). Left: f2°g as a 
function of w^, with the ACDM value marked in green. Right; Effective equation of state lUeff as a function of w^. 



The blue curves in Figure [5] show the dependence of fl'^g and w^g on for e (—3,0). The transition from 
W0 < — 1 to W0 > — 1 is discontinuous since ACDM does not contain dark energy perturbations, and the ACDM value 
is plotted in green. The effective energy density in the constant-w^ models increases monotonically in the quintessence 
regime as u;^ ^ — 1 but asymptotes to fl'^g « 8 x 10^^ in the phantom regime. The effective equation of state acts as 
expected, with that from quintessence being greater than that from ACDM and that from phantom fields being less. 
In both cases the change is relatively slight; for the range £ (—3, —0.15) the effective equation of state is bounded 
by Wee S (0, 1/57) and thus the modifications from a constant equation of state behave approximately as dust. For 
larger equations of state it becomes very negative and acts to accelerate the universe for the rather unrealistic case 
> 18. 

We now turn to quintessence models with a time-varying equation of state. As a first step, we consider two 
examples of basic quintessence potentials, the exponential (see for example [13, [Ml) and the inverse power law [9^ . 
For simplicity and to allow comparison with an EdS case, we do not break the tracking behaviour of the exponential, 
taking 

V^((/.) = Cexp(-a0). (66) 

We choose C such that = 0.8 and = 0.2. Since this field is still tracking the matter component at the present 
epoch, it has an equation of state today of w° w 0. The Ratra-Peebles model we consider instead has a potential 

T/(0) cx 1/02 (67) 

where the energy scale is set such that the energy density has the WMAPV value of il^ = 0.721, which leads to 
w'^ = —0.64. This model accelerates the universe but is illustrative rather than observationally viable. 

Finally, we consider two examples of dynamical dark energy models that satisfy observation constraints. Firstly, 
we employ the standard parameterisation of the equation of state (the 'CPL' parameterisation (93l. [94|). 

w{a) ^ Wo + {1 - a)wa, (68) 

taking wq = —0.9 and Wa = 0.5, consistent with the results in [s^ . Secondly, in order to quantify the influence of 
the scalar field at high redshifts we also consider an early dark energy model. This model is paramcterised in terms 
of the energy density of the quintessence component Qd{a), which evolves in two distinct phases. At low redshifts, it 
behaves as the fractional energy density of a cosmological constant, whereas at high redshifts it contributes a constant 
fraction Q'^ to the total energy density. This can then be written [9^| as 
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FIG. 6: Left: The backreaction in ACDM and selected dynamical dark energy models - the linear CPL parameterisation, a 
Ratra-Peebles field with V{(j)) oc 1/0^, a tracking field with V{(fi) oc exp{—a(t>) and an early dark energy model. Right: The 
effective equations of state for the same models. 

We set the equation of state today wn = —0.99, and choose the fraction of dark energy at early times il^ — 0.02, com- 
patible with observational bounds [96| . The fact that the quintessence component in this parameterisation contributes 
a constant fraction to the total energy density at high redshifts implies an exponential potential for the scalar field 
at those times (with the exponent a in equation (|66p changing as the Universe transitions from radiation to matter 
domination). For low redshifts, the scalar field potential in this model then flattens to produce the desired wq. 



In the left-hand panel of Figure [6] we plot ilcff and AR/R in these dynamical models in comparison with ACDM. 
In contrast with the constant-w^ case, the dynamical models produce an effective energy density lower than that 
of ACDM. The effective energy density at the present day in the CPL parameterisation is ilj^g « 3.8 x 10^*^ and 
the impact on the Raychaudhuri equation is AR/R w —3.6 x 10^^. Backreaction in the model with the Ratra- 
Peebles field is more heavily suppressed with r2|?g « 2.9 x 10~® but the deviation of the Raychaudhuri equation is 
significantly boosted to AR/R « —8.1 x 10~^. The early dark energy model produces H.'^g « 3.3 x 10^^, similar to and 
slightly larger than the Ratra-Peebles form, but generally produces small corrections to the Raychaudhuri equation. 
Finally, the exponential field resembles an EdS model but is heavily suppressed; as we saw earlier, an EdS model with 
these parameters generates ft'^ff « 1.1 x 10^^ at the present epoch, while the exponential potential produces only 

w 7.1 X 10-6. 

The effective equations of state for these models are plotted in the right panel of Figure [S] The CPL, Ratra-Peebles 
and early dark energy models all closely resemble ACDM until a redshift of z w 10 at which point they decline to 
different present-day values, with the Ratra-Peebles field the smallest at w'^g w 1/150, the CPL parameterisation 
dropping to WeS ~ 1/130, and the early dark energy model almost indistinguishable from ACDM. Above z « 10 
the Ratra-Peebles and CPL models are as expected equivalent to ACDM, and the early dark energy model is only 
slightly different. The differences in this case generally arise from Qx>, which is non-zero even in matter domination. 
The effective equation of state produced by the exponential potential naturally displays a very different behaviour 
and closely resembles that from an EdS model suppressed to w^g ~ 1/70. This is because the nature of the dark 
energy perturbations is very different to that of dust perturbations, and as the field is tracking matter it is not diluted 
relative to dust as one moves to higher redshifts. 

Thus far we have only considered dark energy models based on scalar fields with a canonical kinetic term and a sound 
speed of = 1. However, many models of dark energy, such as k-essence, generically predict cl ^ 1 [97l.[98l. [99lll00l |. 
We now investigate the effects of a sound speed different from unity by returning to models with a constant equation 
of state but setting = 0, which models a dark energy allowed to cluster. The dotted lines in Figure[4]show flcs to be 
relatively small for both = —0.9 and = —1.1, while the corrections to the Raychaudhuri equation the deviations 
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are increased further for = —0.9 but suppressed to the level of ACDM for = —1.1.^ The effective equation of 
state for Wfj, = —0.9 is driven back towards the EdS value. The case = —1.1, however, is more interesting. While 
we have allowed the field to "cluster" - although this interpretation is somewhat suspect for a phantom field - the 
effective equation of state is negative, WcS ~ —1/500, although it still closely resembles dust. 

The red curves in Figure [5] show the dependence of and w'^g on when — 0. The dependence of the 
effective energy density on is opposite to the standard case and coincides only at W0 = — 1 where no dark energy 
perturbations are present. For it;^ > —1 the backreaction increases to il'^^ w lO^'^. The "clumps" in the field are 
acting as a dark matter, making this model in this respect resemble an EdS case, although the effect here is somewhat 
lessened. In the phantom regime declines smoothly but increasingly rapidly, again due to the fields dispersing 
structure. The effect is greatly enhanced for = compared to = 1. The effective energy density becomes negative 
at ~ —2.4 and for smaller equations of state becomes increasingly negative. 

In the quintessence regime the equation of state rapidly grows and we recover the EdS value WcS ~ 1/57 only 
for Wfj, « —0.9. For larger values of w^, the effective equation of state increases to a maximum of WcS ~ 1/30 at 
Wcfi = — 1/2 before declining again towards EdS as —^ 0. For the phantom case, however, WcS declines rapidly 
with W0 and for < —1.85, we recover a backreaction which accelerates the universe, albeit with an energy density 
ri[!ff « 10^^ and in a model which is already accelerating. At w —2.4 the effective equation of state passes through 
negative infinity and becomes positive. However, since for these values of the impact on the Friedmann equation 
becomes negative, these models still act to accelerate the universe. 

V. DISCUSSION 

In this paper we have presented spatially averaged cosmological equations in a simple form easily applied to a 
wide variety of models, both numerical and analytic, for which the 3+1 split remains valid and the shift vector is 
vanishing. This encompasses analytic models currently considered in the literature such as LeMaitre-Tolman-Bondi 
models, swiss-cheese generalisations, the Szekeres model or other systems such as networks of McVittie metrics which 
could be used to generate a distribution of sources. The formalism is also applicable to purely numerical approaches, 
the ideal of which would be a fully relativistic n-body code or approximations to such. In our approach, any matter 
source can be expressed as an effective perfect fluid since the physics of the underlying model are left unaltered, and 
so we can consider any combination of fluids ranging from cosmological dust and radiation, to scalar fields, to more 
realistic physical fluids with non-trivial equations of state and couplings between different species. 

Averaging in cosmology has in recent years typically been invoked in an attempt to solve the dark energy problem, 
but is entirely distinct from it and has a history dating back at least to the early 1960s. As such, we have attempted 
to separate the dark energy problem from the averaging problem, and considered systems including not only standard 
matter such as baryons and radiation, along with cold dark matter, but also included scalar fields and dark energies 
on an equal footing. Should scalar fields have played a part in the universe's evolution then it is vital to include them 
in this averaging process. While we have focused on the universe from radiation domination onwards, this statement 
holds true even in the inflationary universe. 

As a concrete example of the impacts one should expect, we applied our formalism to linearly-perturbed Robertson- 
Walker universes. While assuming perturbation theory valid restricts us to minor deviations from Robertson- Walker 
behaviour, the approximation is useful for two chief reasons. Firstly, until recent times the universe is extremely well 
described by linear theory and calculations of the backreaction from linear theory are accurate. Secondly, while the 
results we find for low redshifts may be incomplete they have the advantage of predictive power. A full quantitative 
calculation of the backreaction at recent times is currently unfeasible; as a result we are restricted either to assumptions 
about the nature of the universe (as in [H, H^l), or to perturbation theory (as in [gl, |69|, [7l| ) . Perturbation theory is 
well tested and understood and with modern cosmological probes many of its parameters are known to within a few 
percent. 

We constructed modified versions of the cmbeasy and CMBFast codes to evaluate the backreactions from linear 
modes in a range of models. Far the greatest impacts arise in EdS models. In these models, the effective equation 
of state of the "backreaction fluid" is exactly WcS = 1/57, while the effective energy density is dependent on the 
primordial power spectrum and the Hubble rate at the present day. Changes to the primordial power spectrum 
renormalise the impacts, while changes to the Hubble rate alter them only in matter domination. In radiation 
domination the deviations are independent of the Hubble rate. The effective energy density depends on the Hubble 



^ The coincidence between the dotted curve in Figure|4]and ACDM is, however, just that. Were we to take e.g. = —1.15 we would 
observe a different impact. 
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rate as ^ /i^ ®^^^, and for a standard Hubble rate h = 0.701 we find ft'^g » l O"**. This is approximately 
equivalent to the current energy density in radiation and it has been argued (e.g. jlOl| ) that even small deviations 
can be significant. For a Hubble rate h = 0.45 for which an EdS model can be forced to fit the CMB the impact 
is somewhat less, ilcs ~ 5 x 10~^. The backreaction also evolves in matter domination with the scale factor, and 
so the impact at recombination is of the order of 10^^-10^^, roughly comparable to existing CMB anisotropics. In 
radiation domination the impacts flatten approximately to a constant. The present-day effective energy density for 
the concordance ACDM case is « 4.4 x 10~^, while the effective equation of state undergoes a transition at z « 1 
and diminishes from the pure dust case to Wos ~ 1/120. The effective equation of state from linear modes in these 
two standard cases is thus always greater than zero, in contradiction to BBR08 where the pressure correction St> was 
neglected. The matter-domination prediction of WcS = 1/57 holds until close to recombination. After recombination 
the photon velocities decay and there is a period at which they contribute enough to Tp and 5-p to affect the effective 
equation of state, exhibited as a trough at z « 1000. Before recombination, the effective equation of state increases 
steadily and tends towards WeS = 1/3. 

We then considered a range of dark energy universes, taking as our test cases models with a constant equation of 
state Wff,, the CPL parameterisation u;(a) = wq + {I — a)wa, a quintessence field with an exponential potential, a 
quintessence field with an inverse power law potential, and a model of early dark energy. For the constant equations 
of state we also considered an approximation to a clumping scalar field by setting the field's rest-frame speed of sound 
to zero. 

In the region of = const — 1 with = 1, the dark energy generically creates a slightly larger effective energy 
density and more significant deviations from the usual Raychaudhuri equation than ACDM. The equation of state is 
slightly lower for a quintessence field with > — 1 and slightly higher for a phantom with ly^ < — 1. Allowing 
to vary more widely, the effective energy density tends to vanish as — s- where WeS is marginally negative. For 
realistic models, then, the backreaction remains significantly smaller than the EdS case, and behaves as a dark matter 
with WcB > 0. For strongly phantom cases with Wcf, —3 the effective energy density asymptotes to fi^g « 8 x 10~^ 
and Wcs ~ 1/57. 

When considering dynamical models we set our parameters from WMAPV and as a result selected cases which 
closely resemble ACDM at the background level, the exception being the exponential field which resembles EdS. In 
all cases the backreactions are relatively small. For the CPL parameterisation we recovered backreactions similar in 
both form and magnitude to the ACDM case, although with a slightly smaller present-day effective energy density 
and equation of state. While its effective energy density is very close to ACDM, the inverse power law model produced 
a large correction to the Raychaudhuri equation and an effective equation of state of WcS ~ 1/130. The early dark 
energy model was tuned to closely resemble ACDM at low redshifts, with wq = —0.99 and a dark energy density at 
early times of fi^ = 0.02; unsurprisingly, the backreactions at low redshifts are very close to ACDM although again 
the effective energy density is marginally reduced. At earlier times, the effective equation of state diverges slightly 
from the ACDM case due to a non-vanishing Qx> even in matter domination, although this discrepancy disappears 
before recombination when photon-baryon coupling alters the perturbations significantly. 

Finally, we chose to keep the exponential field as a tracker, to more closely resemble EdS. The backreaction produced 
in this model is surprisingly small, with $1"^ « 7 x 10~^, and the equation of state is reduced to approximately 
Wcff = 1/70 until recombination at which point it again tends towards w^s = 1/3. The reduced effective equation of 
state stems from the different behaviour of the dark energy perturbations; even though the equation of state of the 
field is dustlike, the perturbations differ from dust perturbations. Since the field is tracking matter, fl^/flM does not 
decay for increasing redshift and WcS is thus suppressed at all times. 

In all these cases the effective equation of state is insufficiently negative to accelerate the universe and the back- 
reactions act as dust. This implies that the backreaction on sub-horizon scales acts as a brake on the universe's 
acceleration and a suitably-chosen smaller-scale model might generate a significant impact. This would be similar to 
[ttI where the backreaction from super-horizon modes is used to stop the quintessential expansion. 

As a final case we took fields with a constant equation of state and a sound speed = 0, changing their behaviour 
dramatically. With it;^ = the backreaction is relatively large, fl'^ff ~ 10~^. As this is an approximation of a 
clumping field, the similarity to EdS is perhaps not surprising. Reducing towards Wff, = — 1 rapidly decreases 
the effective energy density of the backreaction. For fields with Wcj, > —1 the effective equation of state reaches a 
peak of WcS ~ 1/30 for « ~l/2, with Wcs ~ 1/57 at « —0.05 and « —0.9. In the phantom regime, the 
effective energy density of the backreaction plummets. For < —1.85, WeS < —1/3 and so the backreaction acts 
to further accelerate the universe. For < —2.4 the effective energy density becomes negative and the equation 
of state positive. Backreaction can then serve to generate a volume-averaged expansion which accelerates, even at 
linear order, albeit at a very low level and in an otherwise unappealing model. It is worth emphasising that while 
the behaviour with = approximates a nonlinear behaviour the magnitude of the effect is a purely linear result. 
Were we to consider a smaller-scale model we might expect a larger backreaction. This may also suggest that for less 
pathological cases with > —1.85, the backreaction could act as a brake on the phantom somewhat. 
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When combined with the results of W03, our study of quintessence cosmologies suggests that further analysis of 
scalar field cosmologies on smaller scales is required. As with any gravitating source, scalar fields are not entirely 
homogeneous. At linear scales we have found that these perturbations generate an effective fluid with Wcs > 
for standard models: re-averaged perturbations in quintessence models do not act to accelerate the universe. If 
we additionally set = to approximate a clumping field, WcS is driven even further from zero, although as the 
calculation is linear the amplitude remains small. The approximate non-linear model of a clumped cosmon in W03 
allowed a better approximation of the amplitude, giving fl^g « 0(1) and w^s ~ —1/15. Since perturbed quintessence 
models generate such high equations of state, and allowing the field to clump alters the behaviour still further, we 
must confirm that an inhomogcneous field behaves on average as a homogeneous field. In forthcoming work we intend 
to tackle this issue employing a direct spatial averaging approach. 

Our approach is not without its issues. The most obvious of these is that due to our methods we have been 
restricted to spatially-fiat universes. This we arc addressing in forthcoming work, again through the use of direct 
spatial averages. More fundamental issues concern both our assumption that the shift vector vanishes and, more 
fundamentally yet, our averaging procedure. The spatial averaging approach averages only scalar projections of the 
Einstein tensor, and these do not contain all the dynamics of the inhomogeneous manifold. It also intrinsically requires 
a 3+1 split, while a more general procedure - Zalaletdinov's macroscopic gravity, for example - can instead perform 
averages in four-dimensional domains. Most directly, we can consider extensions that employ the current formalism 
more generally, the obvious case being nonlinear perturbation theory. While again we would expect only minor 
deviations from standard Robertson- Walker behaviour, this would let us control the scaling of the velocity correctly 
and generalise the nonlinear approximations in BBR08 to non-EdS and non-concordance models. 

In summary, we have presented a general approach to spatial averaging in cosmology applicable to multifiuids in 
a wide variety of metrics, and demonstrated its use with a range of linearly-perturbed Robertson- Walker models. 
This raises interesting questions about the nature of scalar fields, which should be addressed in the near future. We 
have constructed codes which can consider quantitatively any linear Robertson- Walker model. For a model including 
a strongly phantom field with a vanishing speed of sound a backreaction with an equation of state WcS < —1/3 is 
recovered, albeit with a low effective energy density. For standard models the backreaction at linear scales remains 
small and with an equation of state approaching dust. 
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